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A search for the decays B — ► p(770)7 and B° — ► lo(782)j is performed on a sample of 211 million 
T(4S) — ► BB events collected by the BABAR detector at the PEP-II asymmetric-energy e + e~ storage 
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ring. No evidence for the decays is seen. We set the following limits on the individual branching 
fractions B(B + -> p + j) < 1.8 x 10~ 6 , B(B° -> p°j) < 0.4 x 10~ 6 , and B(B° -» wy) < 1.0 x 10" 6 at 
the 90% confidence level (C.L.). We use the quark model to limit the combined branching fraction 
B[B -> (p/w)7] < 1.2 x 10 -6 , from which we determine a constraint on the ratio of CKM matrix 
elements |V t d|/|Vk|. 

PACS numbers: 12.15.Hh, 13.25.Hw 



Within the Standard Model (SM), the decays 
B — > p7 and B° — > W7 proceed primarily through 
a — > c?7 electromagnetic penguin process that con- 
tains a top quark within the loop 0. The rates 
for B+ -> p+7, B° -> p°7, and B° -> W7 are 
related by the spectator-quark model, and we define 
the average branching fraction Q, B[B — > (p/a;)7] = 

i {b(B+ - p+ 7 ) + - A) + B(B° - a^)]}, 

where -^i- is the ratio of B-meson lifetimes. Recent 

T B _ 

calculations of B[B — > (p/a-Oll m the SM indicate a 
range of (0.9 - 1.8) x 1(T 6 0, 0. There may also 
be contributions resulting from physics beyond the 
SM 0. The ratio between the branching fractions 
for B — » (p/uj)j and B — > If* 7 is related in the SM 
to the ratio of Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements |Vtd|/|VjJ 111 Previous searches by 
BABAR 7] and CLEO have found no evidence for 
B — > (p/ui)j decays. 

We search for B — > ,07 and B° — * 0^7 decays in a 
data sample containing 211 ± 2 million T(4S) — > B5 de- 
cays, collected by the BABAR detector :9] at the PEP-II 
asymmetric-energy e + e~ storage ring. The data corre- 
spond to an integrated luminosity of 191 fb _1 . 

The decay B — > pj is reconstructed with p° — > 7r + 7r _ 
and p + — ► 7T + 7r°, while B° — ► luj is reconstructed 
with a; — * 7r + 7r~7r°. Background comes primarily from 
e + e~ — > continuum events, where g = s,c, in 
which a high-energy photon is produced through ir° jr\ — > 
77 decays or via initial-state radiation (ISR). There are 
also significant BB backgrounds: B — > if*7, K* — > Xtt, 
where a if^ is misidentified as a 7T ; B — > (p/uj)tt° and 
B — > (p/u))rj, where a high-energy photon comes from 
the 7T° or 77 decay; and combinatorial background, mostly 
from high multiplicity b — > 57 decays. 

We select 71"^ candidates from tracks with a momentum 
transverse to the beam direction greater than 100 MeV/c. 
The 7r~^~ selection algorithm combines measurements of 
energy loss in the tracking system with any associ- 
ated Cherenkov photons measured by the ring imaging 
Cherenkov detector. The algorithm is optimized to re- 
duce backgrounds from produced in b — > S7 processes 


Neutral pion candidates are identified as pairs of neu- 
tral energy-deposits reconstructed in the Csl crystal 
calorimeter, each with an energy greater than 50MeV 
in the laboratory frame. For B° — > W7 (B + — ► p + j) de- 
cays, the invariant mass of the pair is required to satisfy 



110 < to 77 < 150 MeV/c 2 (117 < to 77 < 145 MeV/c 2 ). 
To reduce combinatorial background, we require the co- 
sine of the opening angle between the daughter photons 
in the laboratory frame be greater than 0.6; this selection 
retains 98% of n° from signal decays. 

A p° candidate is reconstructed by selecting two iden- 
tified pions that have opposite charge and a common 
vertex. We obtain p + candidates by pairing ir° candi- 
dates with an identified tt + . The to candidates are re- 
constructed by combining a ir° candidate with pairs of 
oppositely charged pion candidates that originate from 
a common vertex; the charged pion pair must be con- 
sistent with originating from the interaction region to 
suppress decays. We select p (ui) candidates with 
an invariant mass satisfying 630 < m W7V < 940 MeV/c 2 
(764 < m, +ri p < 795 MeV/c 2 ). 

The high-energy photon from the signal B decay is 
identified as a neutral energy-deposit in the calorimeter. 
We require that the deposit meets a number of criteria 
designed to eliminate background from charged particles 
and hadronic showers [lj}. We veto photons from 7r°(?7) 
decay by requiring that the invariant mass of the candi- 
date combined with any other photon of laboratory en- 
ergy greater than 30 (250) MeV not to be within the 
range 105 to 155 MeV/c 2 (500 to 590 MeV/c 2 ). 

The photon and p/u candidates are combined to form 
the B-meson candidates. We define AE* = E* B — E£ eam , 
where E* B is the center-of-mass (CM) energy of the B- 
meson candidate and £^ cam is the CM beam energy. The 
AE* distribution of Monte Carlo (MC) simulated signal 
events is centered at zero, with a resolution of about 
50 MeV. We also define the beam-energy-substituted 
mass m E s = \J -E^Lm — V*b > where p* B is the CM mo- 
mentum of the B candidate. Signal MC events peak 
in toes at the mass of the B-meson, tub, with a res- 
olution of 3 MeV/c 2 . The distribution of continuum 
and combinatorial BB background peaks in neither Toes 
nor AB*; the background distributions of B — > K*j, 
B — > (p/u))ir° and B — > (p/u>)r) peak at tub in toes and 
between -190 MeV and -60 MeV in AE* . We consider 
candidates in the ranges —0.3 < AB* < 0.3 GeV and 
5.20 < toes < 5.29 GeV/c 2 to incorporate sidebands 
that allow the combinatorial background yields to be ex- 
tracted from a fit to the data. 

Several variables that distinguish between signal and 
continuum events are combined in a neural network |llj | . 
The input variables depend mainly on the rest of the 
event (ROE), defined to be all charged tracks and neu- 
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tral energy deposits in the calorimeter not used to re- 
construct the B candidate. To reject ISR events, we 
compute the ratio of second-to-zeroth order Fox- Wolfram 
moments f° r the ROE and the p/uj candidate, in 
the frame recoiling against the photon momentum. To 
discriminate between the jet-like continuum background 
and the more spherically-symmetric signal events, we 
compute the angle between the photon and the thrust 
axis of the ROE in the CM frame and the moments 
Li = J2jP*j ' I cos6j\ % / J2jPji where p* and 9* are the 
momentum and angle with respect to an axis, respec- 
tively, for each particle j in the ROE. We use L\, L2, 
and I/3 with respect to the thrust axis of the ROE, as 
well as L\ with respect to the photon direction. Differ- 
ences in lepton and kaon production between background 
and B decays are exploited by including BABAR flavor 
tagging variables 01 as weu as the maximum CM mo- 
mentum and number of and K® in the ROE. For the 
B° — > (p°/uj)j modes, we also use the separation along 
the beam axis of the B-meson candidate and ROE ver- 
tices; to remove poorly reconstructed events we require 
the separation be less than 4 mm. A separate neural 
network is trained for each mode. We make a loose selec- 
tion on the output of the neural network, J\f, that retains 
around 80% of the signal events. 

To suppress background, we combine a number of 
signal-decay variables in a Fisher discriminant [l^j . 
separately for each mode. We calculate the _B-meson 
production angle 9* B , the p/cu helicity angle 8h, which is 
defined with respect to the normal of the decay plane 
for u) candidates, and the u> Dalitz angle 9d |7j. To 
reject B — > p(TT°/r]) and B — > uj(tt° /rj) events in the 
B + — ► p + 7 and B° — ► uj^ (B° — * p°"f) selection, we 
require | cos9 H | < 0.70 (0.75). 

After applying the AT and | cos Or | criteria the expected 
average candidate multiplicity in signal events is 1.15, 
1.03 and 1.14 for B+ p+j, B° -> p°~f and B° -> ujj, 
respectively; in events with multiple candidates the one 
with the smallest value of |A_E*| is retained. 

The signal yield is determined from an extended max- 
imum likelihood fit to the selected data. We fit the 
four-dimensional distribution of toesi AE* , T and Af. 
For the B — > pj fits, five event hypotheses are con- 
sidered: signal, continuum background, combinatorial 
-B-background, peaking B — > p{-K° /rj) background and 
peaking B — > K*j background. For the B° — > ui-f fit we 
consider only signal, continuum background, and peak- 
ing B — > u)(ir°/r]) background. The correlations among 
the observables are small; therefore, we assume that the 
probability density function (PDF) V{xj\ ai) for each hy- 
pothesis is the product of individual PDFs for the vari- 
ables Xj = {mES, AE* , T, A/"} given the set of parameters 
d?j . The likelihood function is a product over all Nk can- 



didate events of the sum of the PDFs, 



C k = exp 



E 



n 



'Ni 



hyp 



ri'il^i (Xj 



where rij is the yield of each hypothesis, k is the B — > 
(p/w)7 mode, and Ah yp = 5(3) for B — > fry (B — ► ory). 

The ttt-es and AE* PDFs are parameterized by a Crys- 
tal Ball function 0] for both the signal and peaking 
background. The parameterization is determined from 
signal MC samples, except the mean of the Ai?* dis- 
tribution, which is offset by the observed difference be- 
tween data and MC samples of B — > K*-f decays. The 
continuum background wies and AE* distributions are 
parameterized by an ARGUS threshold function and 
a second-order polynomial, respectively. The combinato- 
rial B background is described by a smoothed distribu- 
tion [l7^ determined from MC events in both ibes and 
AE*. The distribution of Af for signal and BB back- 
ground is parameterized by a Crystal Ball function. The 
Af distribution for continuum is determined from side- 
band data, and a histogram is used as the PDF. The dis- 
tribution of J- is parameterized by smoothed histograms 
of sideband data for the continuum background and MC 
events for all other hypotheses. 

The fit to the data determines the shape parameters 
of the continuum background toes and AE* PDFs, as 
well as the signal, continuum background and combina- 
torial BB background yields. All other parameters are 
fixed, including the peaking BB background yields. A 
combined fit is also performed relating the modes using 
the definition of B[B — > (p/uj)j} to determine an effective 
yield (n e s) assuming n(B + — > p +r y) = n c fi-e(B + — > p +r y) 
and n(B° -> (p%>) 7 ) = l^eff • <B° - (jP/wfr) 
where n and e are the yields and reconstruction efficien- 
cies of each mode; the efficiencies include the daughter 
branching fractions. We take r -f^ = 1.086 ± 0.017 0. 
Fig. H shows the projections of the combined fit results 
compared to the data. The results for the individual 
mode signal yields and n s are given in Table [I] The 
significance is computed as \J1 A log £ where Alog£ is 
the log likelihood difference between the best fit and the 
null-signal hypothesis. No significant signal is observed. 

The most important systematic uncertainties are asso- 
ciated with the modeling of BB backgrounds, the fixed 
parameters of the PDFs used in the fit, and the signal 
reconstruction efficiency. The first two contribute to the 
uncertainties on the signal yields. 

The uncertainty on the peaking B — > K*-f back- 
ground is dominated by the misidcntification rate; 
the rate is corrected by the difference in AT* misiden- 
tification between data and MC samples of D* decays, 
with the whole correction taken as the uncertainty. For 
the B+ -> p+(7r%), B° -► p°r) and B° -> c^tt /^) 
peaking background decays, we vary the branching frac- 
tions by either one standard deviation from the measured 
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TABLE I: The signal yield (n s i g ), continuum background yield (n cont ), peaking background (n pea k), significance in standard 
deviations a, efficiency (e), and branching fraction (B) central value and upper limit at the 90% C.L for each mode. The results 
of the combined fit are shown in the bottom row where n s i g is equal to n e ff , which is described in the text. When two errors 
are quoted, the first is statistical and the second is systematic. 



Mode 


Uslg 




Significance (a) 


e(%) 


S(10" 6 ) 


B(HT 6 ) 90% C.L. 


B+ -> p+7 


nc+15+2 
zo -14 -2 


6850 ±90 18 ±4 


1.9 


13.2 ± 1.4 


0.9 t 0.5 ± 0.1 


< 1.8 


B° - A 


n Q + 7.2 + 1.7 

u -°-5.4-1.6 


4269 ±73 18 ±7 


0.0 


15.8 ± 1.9 


0.0 ±0.2 ±0.1 


< 0.4 


B° -s-ary 


o q+5.7+1.3 
°'°-4.5— 1.8 


1378 ± 37 2.6±?;| 


1.5 


8.6 ±0.9 


0.5 ±0.3 ±0.1 


< 1.0 


Combined 


9RQ 4-l2ti+40 
ZU»_ 12 Q_4 5 




2.1 




0.6 ±0.3 ±0.1 


< 1.2 




1.2 5.22 5.24 5.26 5.28 
m ES (GeV/c ) 




1.3-0.2-0.1 0.1 0.2 0.3 

AE* (GeV) 




0.7 0.8 0.9 1 

J\f 



FIG. 1: Projections of the combined fit to B — ► pj and B° — + 
uj'y in the four discriminating variables: (a) weSi (b) AE* , 
(c) TV, and (d) T . The points are data, the solid line is the 
total PDF and the dashed line is the background only PDF. 
The selections applied, unless the variable is projected, are: 
5.272 < m ES < 5.286 GeV/c 2 , -0.10 < AE* < 0.05 GeV 
and M > 0.9; the selection efficiencies for signal events are 
45%, 57%, 70%, and 44% for the m ES , AE* , Af and T pro- 
jections, respectively. 



values or between zero and the measured upper limit if 
the decay has not been observed 0, |2(| ; the value of the 
B° — > p°ir° branching fraction is varied between zero and 
5.1 x 10~ 6 H(l El . The uncertainty on the peaking back- 
ground of each mode is shown in TableQ] We find the bias 
from neglecting the B — > K*j background and combina- 
torial BB background in the fit to B° — > ujj candidates 
is 1.1 jli'i events; the corrected yield is given in Table [Q 
To estimate the uncertainty related to the extraction of 
the signal tties and AE* PDFs from MC distributions, 
we vary the parameters within their errors. The variation 
in the fitted signal yield is taken as a systematic uncer- 
tainty. The uncertainty related to the statistics of the 
histogram PDF that describes the continuum J\f distri- 
bution is evaluated by varying the binning and by using a 



fifth-order polynomial as an alternative PDF. Several dif- 
ferent control samples of data and MC events were used 
to determine alternative PDFs for the different hypothe- 
ses; none of these resulted in a significant change to the 
fitted signal yield. 

The signal efficiency systematic error contains uncer- 
tainties from tracking, particle identification, photon/7r° 
reconstruction, photon selection and the neural network 
selection that are determined as in Ref. j^. We deter- 
mine the effect of correlations among the fit variables by 
using an ensemble of MC experiments of parameterized 
continuum background simulations embedded in samples 
of fully simulated signal and BB background events. No 
bias is observed within the statistical error on the mean 
yields from this ensemble, which is taken as a multi- 
plicative systematic uncertainty. The total multiplica- 
tive systematic error values are 11%, 13% and 10% for 
B + — ► p + 7, B° — ► / o°7 and B° — > u>j, respectively. The 
corrected signal efficiencies and their uncertainties are 
shown in Table [I] 

In calculating branching fractions, we assume 
B{T{AS) -» B°B°) = B(T(4S) -> B+B~) = 0.5. The 
90% C.L. is taken as the largest value of the efficiency- 
corrected signal yield at which 2Alog£ = 1.28 2 . We 
include systematic uncertainties by increasing the effi- 
ciency corrected signal yield by 1.28 times its systematic 
uncertainty. Tabled shows the resulting upper limits on 
the branching fractions. 

Using the measured value of B(B — > K*j) j^, we cal- 
culate a limit of B[B -> (p/u)~/]/B{B -> K*j) < 0.029 
at 90% C.L. This limit is used to constrain the ratio_pf 
CKM elements |Vtd/Vts| by means of the equation 



B[B -> (p/uh) 
B(B -> K*-/) 



V t 



id 



1 - mg/Ag 
,/Mj 



1 



C 2 [l±Ai?], 



K 



where C describes the flavor-SU(3) breaking between p/uj 
and K*, and AR accounts for annihilation diagrams. 
Both £ and AR must be taken from theory 0, 0, l23| . 
Following 0], we choose the values C = 0-85 ± 0.10, 
and AR = 0.10 ± 0.10, which is the average over the 
values given for the three modes. We find the limit 
|V trf |/|y ts | < 0.19 at 90% C.L, ignoring the theoreti- 
cal uncertainties. Our upper limit on |Vtd|/|Vt a | con- 
strains iVtdl < 0.008 at 90% C.L. assuming \V ts \ = \V cb \ 
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[lij : this lies within the current 90% confidence interval 
0.005 < |Vtd| < 0.014, which is obtained from a fit to 
experimental results on the CKM matrix elements 
Varying the values of £ and AR within their uncertain- 
ties leads to changes in the limits by ±0.03 and ±0.001 
for IVfdl/IVfal and \V td \, respectively. 

In conclusion, we have found no evidence for the ex- 
clusive b — > c?7 transitions B — > pj and B° — ► oj~f in 211 
million T(4S) -> SB decays studied with the BABAR 
detector. The 90% C.L. upper limits on the branching 
fractions and |Vtrfl/|Vt a | are significantly lower than our 
previous values Jbland restrict the range indicated by 
SM predictions UH- 
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